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We calculate the conductivity of a clean graphene sheet at finite temperatures starting 
from the tight-binding model. We obtain a finite value for the dc-conductivity at zero 
temperature. For finite temperature, the spontaneous electron-hole creation, responsible 
for the finite conductivity at zero temperature, is washed out and the dc-conductivity 
yields zero. Our results are in agreement with calculations based on the field-theoretical 
' O ' model for graphene. 
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1. Introduction 



The recent achievements of isolating and locating a single layer of graphene and 
performing transport measurements have tremendously stimulated the research on 
l/^ . carbon-based electronics.^^ By applying an external gate voltage, the system can be 

switched from n-type to p-type carriers, thus controlling the electronic properties 
and giving rise to carbon-based devices. At the crossover from n- to p-type carri- 
, ers, i.e., at the charge neutrality point where no carriers are present, experiments 

nevertheless report a minimal finite, "universal" conductivity.!^ 

A first understanding of this result is already obtained from the band structure of 
the system. Within a two-band tight-binding approacfPI and in the long- wavelength 
limit, single-layer graphene is described by two-dimensional Dirac Fermions (see 
. e.g. Ref. At zero doping, the density-of-states (DOS) at the Fermi energy is 

zero and the system is a zero-band gap semiconductor or semi-metal. Graphene 
thus lies at the borderline of a semiconductor and a metal which gives rise to a 
number of new phenomena. E.g., in contrast to the zero (infinite) conductivity for 
a clean semiconductor (metal), the conductivity of a clean graphene sheet at zero 
temperature is finite and of order e'^/h. 

Already Ludwig et. al.l^l pointed out that there are two possible prefactors to the 
conductance quantum e'^/h. Calculating, e.g., the dc-conductivity via the Kubo- 
formula using the polarizability for undoped graphene of Ref. ^ yields = fx' 
where the dielectric constant e is usually omitted (e — 1). Introducing a finite 
scattering amplitude F due to impurities and performing the limit F — > in the end 
leads to a ^^4^ 
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Besides using the Kubo formula, there were also calculations based on the Lan- 
dauer formula for a rectangular system with aspect ratio W/L ^ l P l m i Independent 
from the explicit boundary conditions, the conductivity was given hj a = cor- 
responding to the physical limit procedure (first w and then F ^ 0) in the 
case of the infinite system without contacts. In Ref. ^I^^, Ryu et. al. thoroughly dis- 
cuss the crossover between the two values depending on the different dc-limits. But 
taking the dc-limit after the integration over energies is performed, one obtains a 

2 

third finite value for the minimal conductivity, i.e., a = tt^ which is closest to the 
experimentally measured value given by ct « 4^.^!^ 

The above works are based on a ballistic transport model which predicts a uni- 
versal conductivity. But Mishchenko pointed out that electron-electron interaction 
may lead to a non-universal correction.H^I Also diffusive transport m odel s can ac- 
count for a minimal, but non-univers al c onductivity via midgap states,!^ Coulomb 
scattererJlS! or electron-hole puddles .^1^ 

This paper shall focus on the transport properties of clean graphene sheets at 
finite te mp erature and frequency. For an overview on ac conductivity in graphene 
see Ref. UJ. The effect of temperature was also discussed in Ref. which can have 
two contrary effects as outlined below. 

On the one hand, one can see temperature as some effectiv e fi nite chemical po- 
tential which would lead to an enhancement of the conductivity.^^ This is suggested 
by the results by Vafek on the plasmon dispersion for graphene at zero doping and 
at finite temperature He finds that the collective plasmon excitations are only 
weakly damped even though they lie in the region where electron-hole excitations 
are allowed. In the case of finite doping and zero temperature, electron-hole excita- 
tions are forbidden due to the Pauli principle, thus leading to undamped plasmon 
excitations within the RPA approximation.'^'-' ^ll 

On the other hand, the finite conductivity at T = can be related to the 
spontaneous electron-hole creation at zero energy.I^l These quantum fluctuations 
can be washed out at finite temperature, thus rendering the mechanism for finite 
conductivity for clean graphene to zero. It is this scenario which prevails and the 
dc-conductivity for graphene is indeed zero if the energy scale set by the carrier 
density is smaller than the temperature. 

All calculations on the conductivity of clean graphene cited above are based 
on the long- wavelength approximation, i.e., start from relativistic Dirac Fermions. 
Here, we want to discuss the conductivity of a clean graphene sheet starting from 
the tight-binding model. We will further include finite temperature since the order 
of limits Lo and then T — > corresponds to the physical realization of an 
experiment. 
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2. Tight-binding model with vector potential 

The simplest Hamiltonian describing non-interactive electrons on graphene in the 
presence of a time dependent vector potential A{t) reads 

H = -J2 [tR,R+6al{R)b,{R + S) + H.c] , (1) 

where R runs over all unit cells and S runs over all three nearest neighbors with 

51 = ^{-l,V3,0), (2) 

52 = |(-l,-y3,0), (3) 

53 = a(l,0,0). (4) 

The bare tunnel-matrix element t is thus modified by the time dependent vector 
potential A{t) as 



/ e f^+^ \ / e \ 

tR,R+s = iexp I «27r- J dl ■ A{t) 1 = texp \J,2Tr -A{t) ■ Sj . 



(5) 



The current operator in the presence of the vector potential is defined as 

j = -VAH. (6) 

In linear response one needs the current operator up to first order in A[t), which 
gives 

j = tYl ( "^^^ + Uydy)al{R)ba{R + S) 

R,a,d 

( u,AJl + UyAySl)aU-R)b^i-R + 6)+ H.c.) . (7) 

The Maxwell equation 

allows for a simple relation between the vector potential and the electric field, 
which is needed for the calculation of the electrical conductivity. If one introduces 
the Fourier representation for the operators as 

a<.(i«) = ^^e-'=-«a<.(fe), (9) 



and the equivalent expression for the operators b{R), the Hamiltonian, for A = 0, 
reads 

H = -tJ2[Hk)al{k)b4k) + R.c.], (10) 

k,a 
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where 

s 

It proves useful to redefine (j){k) as 

(/.(fc) = 5]e-*'=-(''-''^), (12) 



and introduce the transformation After these transformations, the 



4 ^ e-'''-^'al. 

single particle Green's functions for the Hamiltonian are given by 

^'Asi^n, k) 

where we also introduced the chemical potential fj,. 
3. Current operator and Kubo formula 

Let us now concentrate on the current operator jx, which is composed of the para- 
magnetic and diamagnetic contribution — + Ax{t)j^, each of them given 

by 

= Y.^im - i)ai{k)K{k) - (r (fc) - 3)6t , (17) 





iuJn + M/^ 


(13) 








-t(j)*{k)/h 


(14) 


(iuJn - 


^fi/hy^t^\4>{k)\yh^ ' 






(15) 


{iuJn - 


^fi/hy-t^\ci){k)\yh^ ' 




-t(l){k)/h 


(16) 







2h 

k., 



and 



jf - T^^im + i)al{k)K{k) + (0*(fe) + 3)6t (fe)a,(fe)] . (18) 

The Kubo formula for the conductivity is given bjI^H 

I ^_ <jf > , Ax.(^ + »0+) .^„. 
''""^'^^ " iA,{uo + iO+) ^ ihA,{uj + iO+) ' ^ ^ 

with As = NcAc the area of the sample, and Ac the area of the unit cell, from which 
it follows that 

where D is the charge stiffness which reads 

hA. • (2^) 
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The incoherent contribution to the conductivity h.^xi'^ + iO^) is obtained from 
Kxx{i^n), with this latter quantity defined as 

.hp 

kxxiiu^n)^ dre'-"-<T.jf(T)jf(0)> . (22) 







The /—sum rule, giving the oscillator strength, reads 







axxiuj)duj = -TT^^ — , (23) 



where we have used 6{u!)dLu = 1/2, see also Ref. 1221, 
Within this model, the several quantities read 

< E \m\[nF{-t\m\ - - nF{t\m\ - m)] , (24) 

k 

^xxiO + ^0+) = ^ E ll^f ["^(-^l'^^^)! - ^) - nprnk)] ^,)] , (25) 
QAxxico + iO+) = J2 f[m][nF{-t\m\ - - nFit\^ik)\ - /.)] 

k 

X [nS{uj - 2t\(j){k)\/h) - nS{uj + 2t\(j){k)\/h)] , (26) 

and 

/l«.,|.i8-4w.,p + i8MMy|*ffl!!^ p„ 

The above formulas were derived using the fact that 



\j2\^ik)fgi\m\) (28) 



3 

k 

where g{\(j){k)\) is an arbitrary function of the absolute value of 0(fc), and the fact 
that in the Dirac cone approximation one has 

1^ ^ e^'-/' [cos(20) - ^ sin(20)] , (29) 

and a similar expression for the complex conjugate expression [(p* (k)]'^ /\(j){k)\'^ . One 
can also show that the following relation holds true 

which proves that the charge stiffness is zero at zero temperature and for zero 
chemical potential (half-filling). As a consequence, the system can only show d.c. 
conductivity at half filling and zero temperature if 

CTj „ = hm — (31) 
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is finite. 

The calculation is simple to do within the Dirac cone approximation. The term 
in f[<f>{k)] proportional to — [50(fc)]^ gives zero in this approximation. Let 

us first introduce the density of states per unit area p{e) 



1 

4^ 



qc 3 

27rqdqS{e — -aq) 



4e 



187ra^ 



(32) 



Using p{e) we can calculate 
1 



A. 



■ Y,f[m]sii^-2t\m\/fi) 



dep{e){18 - 4:e^)6{uj - 2te/h) 



(33) 



The contribution just computed corresponds to the value of a single Dirac cone. 
The optical conductivity is (the two Dirac cone contributions included) 



cr(w) 

which gives a d.c. value of 



TT e 
2~h 



\3V2t) 



(34) 



vr e 
2~h 



(35) 



This is the result which one obtains if a finite damping term is not included.^m 

The above result holds only at zero temperature. For finite temperature, which is 
the case in any experiment, the conductivity, at half filling and within the Dirac cone 
approximation, is (note that due to particle-hole symmetry the chemical potential 
is zero for any temperature) 



leading to 



^ (y3V2t) 



ad.c.(r) = o, 



tanh 



Wbt) ' 



(36) 



(37) 



which should be interpreted as the correct result at zero frequency instead of result 
(|35p . This is shown in Fig.[Tl where the optical conductivity as function of frequency 
is plotted for various temperatures for t — 3eV. Only the T = curve yields a finite 
dc conductivity whereas for finite temperature all curves eventually become zero 
for u! 0. E.g. at T = 70K, the universal value is only reached for frequencies 
uj > 0.075eV, at T = 4ii: for uj > O.OleV. 
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Fig. 1. The optical conductivity as function of frequency is plotted for various temperatures (for 
t = 3eV). 



4. Summary 

Starting from the tight-binding model, we have calculated the conductivity of clean 
gaphene as function of the frequency and temperature. Depending on the order 

2 

of the limits uj ^ and T — > 0, we obtain either the finite value o" fx 
or zero in agreement with field-theoretical calculations \^ The universal value of 
the conductivi ty o f a clean graphene at finite frequencies was recently confirmed 
experimentally.'^ 

Our calculations start from the tight-binding model to complement the current 
discussion on the conductivity of clean graphene and we believe that the explicit 
calculation based on the lattice model will be useful to a wide audience. An obvious 
extension is to include the next-nearest neighbor coupling and corrections to the 
Dirac cone approximation.^^SI Also, the effect of out-of-plane phononJ^Sl ^nd other 
scattering mechanismJ23 on the optical conductivity of clean graphene can easily 
be assessed within the presented formalism. 
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